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This paper is an English translation of Michel Henon's article, Sur 
revolution dynamique des amas globulaires. II- L'amas isole originally pub- 
lished in French in the Annales d'Astrophysique, Vol. 28, p. 62 (1965). 

A translation of the first paper of this series (Paper I: Sur revolution dy- 
namique des amas globulaires, Henon 1961, Annales d'Astrophysique, Vol. 24, 
p. 369) is also available. 

Conventions and notations are as in the original version, for consistency. 
The English version is written so that it is as faithful to the French text as 
possible. The translator added some notes [T.N] for the sake of clarity, when 
required. French, English and Russian abstracts written by M. Henon are 
available in the original version of the paper. The English part is reproduced 
below. 

FR thanks Michel Henon for the enthusiasm and kindness he expressed 
when he was asked for permission to translate his work, as well as Douglas 
Heggie and Mark Gieles for their careful proofreading. 
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ORIGINAL ABSTRACT 



The homologic [T.N.: homologous] mod el which was 



defined and studied in an earlier paper (jHenonI 119611 ) 
is computed for the case of an isolated cluster. This 
model is compared to von Hoerner's numerical results 
( von Hoernerll963h . The agreement is quite good for the 
density low [sic law]; the distribution of energies and the 
rate of escape. The broad features of the evolution are: 
tendency towards the homologic model, building up of an 
infinite central density, formation of close binaries, expan- 
sion of the cluster as a whole. An unexplained discrepancy 
remains in the rate of evolution. 

1 INTRODUCTION 



In a previous work (|Henonlll961l . hereafter Paper I), we 
have obtained a system of fundamental equations that 
allows one to calculate the dynamical evolution of a star 
cluster from a given initial state. We have then calculated 
and studied in detail a particular solution of this system: 
the "homologous model" , which has the property of re- 



maining self-similar in time, the evolution being reduced 
to the scaling of the physical quantities. 

The cluster was supposed to be plunged into an inho- 
mogeneous field of forces, due to the Galaxy. The effect of 
such a field (often called "tidal effect" ) is to limit the clus- 
ter radius to a specific value r^, linked with the mass of 
the cluster through the relation (Paper I, Equation 3.8): 



re oc 



(1) 



Every star further than the distance re definitely es- 
capes from the cluster. This way, a loss of stars is con- 
tinuously happening, and the cluster mass decreases with 
time. From Equation 1, the radius decreases too. 

The present paper focusses on the different case of an 
isolated cluster^ i.e. when the external forces are zero or 
negligible. Then, the radius Tg is infinite, and the stars 
cannot escape anymor43 (see Paper I, Chapter HI). Here- 
after, the mass A^e of the cluster is constant over time. 
All the equations given in Paper I are still valid in this 
case, as soon as A = cx), and one can solve them by use of 
to numerical computations (see Paper I, Chapter V). 



more detailed study of this phenomenon is given inlHenorJ l| l96Cj) ; we find that the escape rate is not exactly zero, but negligible in 
practice. Similar results have been obtained by WooUcv fc Dickcng ||1962| ). 
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Figure 1: Isolated homologous model: the four fundamental func- 
tions F(E), the distribution function; D(U), the spatial density; 
R(U), the distance to the centre; Q(E), the adiabatic invariant. 
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Figure 2: spatial density of the isolated homologous model (on top 
and right) and of the non-isolated model (on bottom and left), as a 
function of the radius. 



2 RESULTS 
2.1 Structure 

Once again, we find that only one unique homologous 
model exists. The four fundamental functions F(E), 
D(U), R(U), Q(E) are given in Table [2 and plotted in 
Figure [TJ As in Paper I, all the variables relative to the 
homologous model are in bold font; E is the energy, U is 
the potential, F is the distribution function, D the spatial 
density, R the distance to the centre; Q is a functional 
form without particular physical meaning (one can show, 
in fact, that Q is the adiabatic invariant associated to the 
total energy E, i.e. the volume of phase-space en closed 

Fig- 
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within the surface E — est; see 
ure[T]is quite similar to the analogous figure in Paper I, 
the main difference being that here, R goes to infinity 
when U goes to zero, because the cluster has an infinite 
radius. 

Below E — — 5orU — —5, one can use the asymptotic 
forms of (4.34) and (4.35) from Paper I. Furthermore, for 
E — >■ 0, or U — > 0, the asymptotic forms are: 



F 
D 
R 

Q 



flo e 



(2) 



= 2V^, 



-3/2^_U)9/'* 



-Me/U, 

^^(-E) 
12 V2 



-3/2 



with: 

ao = 6.391, a = 3.544. 
The values of constants found in the equations are: 



K 

c 

b 

Me 



1.363, 
-1.097, 



(3) 



1.946 
1.217 



-0.7313, 

(total mass), 

(total kinetic energy). 



The spatial density D is plotted in Figure [5] as a func- 
tion of the radius, and compared to those of the non- 
isolated model (the two curves are normalized to get the 
same asymptotical form for R — > 0). Toward the out- 
skirts, the models are naturally very different, as one 
yields a finite radius while the other is infinite. How- 
ever, close to the centre, both models overlap. Therefore, 
the presence of an external field only modifies the outer 
layers of the cluster, as one could have expected. To be 
clearer, we can state that the two models are almost iden- 
tical within the sphere enclosing half of the total mass 
(the radius of this sphere, the so-called median {[T.N.: 
half-mass]) radius, is R = 0.7563 for the isolated cluster). 
In particular, we note that the isolated cluster, like the 
non-isolated cluster, yields a density that goes to infinity 
as R~^ when R goes to zero. 

The last column of Table[T]gives the projected density 
Dp as a function of R. It goes to infinity as R~^ when 
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Table 1: The isolated homologous model 



E or U 


F 


D 


R 


Q 


Dp 


-5 


131.81 


134.66 


0.1126 


0.0002927 


39.86 


-4.8 


106.43 


106.45 


0.1256 


4025 


34.72 


-4.6 


85.80 


83.83 


0.1402 


5554 


30.13 


-4.4 


69.05 


65.71 


0.1566 


7680 


26.05 


-4.2 


55.46 


51.26 


0.1753 


0.001064 


22.41 


-4 


44.44 


39.75 


0.1964 


1479 


19.18 


-3.8 


35.52 


30.63 


0.2204 


2060 


16.33 


-3.6 


28.30 


23.43 


0.2478 


2879 


13.80 


-3.4 


22.47 


17.76 


0.2793 


4039 


11.58 


-3.2 


17.77 


13.34 


0.3153 


5691 


9.629 


-3 


13.98 


9.898 


0.3571 


8058 


7.928 


-2.8 


10.94 


7.247 


0.4058 


0.01148 


6.451 


-2.6 


8.495 


5.220 


0.4629 


1647 


5.176 


-2.4 


6.544 


3.689 


0.5304 


2382 


4.085 


-2.2 


4.989 


2.546 


0.6110 


3479 


3.160 


-2 


3.754 


1.707 


0.7085 


5141 


2.384 


-1.8 


2.777 


1.104 


0.8283 


7709 


1.744 


-1.6 


2.010 


0.6818 


0.9782 


0.1177 


1.226 


-1.4 


1.413 


0.3958 


1.171 


1837 


0.8170 


-1.2 


0.9516 


0.2112 


1.425 


2955 


0.5060 


-1 


0.6018 


0.09978 


1.778 


4948 


0.2816 


-0.9 


0.4623 


0.06429 


2.011 


6533 




-0.8 


0.3437 


0.03904 


2.300 


8779 


0.1324 


-0.7 


0.2444 


0.02191 


2.667 


1.206 




-0.6 


0.1631 


0.01103 


3.151 


1.706 


0.0462 


-0.5 


0.09898 


0.004740 


3.822 


2.510 


0.0227 


-0.4 


0.05150 


0.001584 


4.817 


3.909 


0.00928 


-0.3 


0.02015 


0.000350 


6.457 


6.66 


0.00262 


-0.2 


0.00371 


0.000033 


9.710 


13.4: 


0.000346 


-0.1 


0.00011 


0.000000 


19.41 


41.1: 


0.000005 
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R goes to zero. For R — > oo, its asymptotic form is: 

Dp = 4^0 <J-^ R-3/2g-.VR7^. (4) 

2.2 Evolution 

The equations (5.14) of Paper I, which give the two pa- 
rameters /3 and 7 of the homology as functions of time, 
become: 



1 



70 



T 

1 



-2/3 



(5) 



From this, one can derive the variations of the physical 
quantities as functions of time. In particular, the radius 
of the cluster evolves as 



2/3 



(6) 



Here, r represents any quantity which characterizes the 
size of the cluster, e.g. the median radius defined above. 
ti is negative (see Paper I, 5.13). The variation of r is 
plotted on Figure [3] We see that the isolated cluster 
expands with time (while the non-isolated cluster con- 
tracts). 



n. - 




Figure 3: Isolated homologous model: evolution of the radius. 



This result differs from the usual idea according to 
which internal encounters make the cluster contract. This 
can be easily explained. The classical reasoning is as fol- 
lows: the stars which escape from the cluster take with 
them practically zero energy, so that the total energy of 
the cluster remains constant; according to the virial the- 
orem, the potential energy also remains constant; as the 
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total number of stars decreases, the distances between 
them must decrease too. But here, on the one hand there 
are no escapers because the cluster is isolated. On the 
other hand, the fact that the parameter K is non-zero 
and negative (as shown in Paper I, Chapter IV) implies 
the existence of a negative flux of energy toward the cen- 
tre of the cluster. The phenomenon occurs in the same 
way as for non-isolated clusters. This negative energy 
accumulates in the centre, thanks to the formation of in- 
creasingly close binaries or multiple stars. At the same 
time, the energy in the rest of the cluster decreases (in 
absolute value); inducing an increase of the distance be- 
tween the stars, i.e. a dilation of the cluster. 

The curve on Figure [3] continues to infinity in the fu- 
ture; it is therefore impossible to define a remaining life- 
time for the cluster, as opposed to the non-isolated case. 
In theory, the evolution of the cluster goes on indefinitely. 
In practice however, while the cluster expands continu- 
ously, it will eventually feel the gravitational effect of an- 
other object; we are back to the non-isolated case studied 
in Paper I. 

The curve in Figure [3] is limited toward small i's. 
Thus, one can obtain an upper limit equal to |ti| for the 
age of the isolated cluster. This is strictly true only when 
one assumes that the cluster has always followed the ho- 
mologous model, but is probably true in an approximate 
fashion in other cases. This limit can be calculated for a 
real cluster of which one knows the mass and size, thanks 
to (2.24), (5.13) and (8.18) of Paper I; one obtains: 



1 



Ql/2 



In (n) 



1/2 



3/2 



(7) 



3 COMPARISONS WITH THE NUMERICAL 
EXPERIMENTS OF VON HOERNER 

There is no such thing as a real system with which the 
present model could be usefully compared. The open clus- 
ters and the globular clusters can surely not be considered 
as isolated; their observed radius is close to the limit set 
by the galactic field, in general. The elliptical galaxies 
could, in many cases, be considered as isolated systems 
to a good approximation; but their internal relaxation 
time is extremely long, so that their structure is probably 
not set by [T.N.: stellar] encounter effects. Finally, for 
galaxy clusters, the distribution function cannot be pre- 
cisely determined, in particular in the external regions. 

However, the present model ca n be compare d with 
the two "experimental clusters" of von Hoerner (|l963l 
hereafter v.H.). These are imaginary clusters made of 
25 stars, whose evolution has been numerically computed 
for a time much longer than the relaxation time. No 
other forces than mutual attraction are taken into ac- 
count; therefore, these clusters lie in the isolated case. 
Furthermore, the masses of the stars are taken to be 
equal; as in the homologous model. 



3.1 Density law 

The evolution of the spatial distribution of the clusters is 
plotted in v.H., Figure 5. For the comparison, we con- 
sider the final state, the most relaxed, i.e. t — 5.4 in the 
first {[T.N.: experimental cluster]) case and t — 2.9 in 
the second case (corresponding to the mean of the fifth 
and fourth last columns of Table 2 of v.H., respectively). 
In order to compare these distributions with those of the 
homologous model, it is necessary to scale r and p cor- 
rectly. These scalings derive from the median radius and 
the total mass. The two clusters have median radii of 
0.840 and 1.006 respectively (from v.H., Table 10) and a 
total mass of unity; the homologous model has a median 
radius of 0.7563 and a mass of: 

/"OO 

/ 47rD R2 dR = 47rMe = 24.45 (Paper I, 2.45). 
Jo 

(8) 



One finds: 



R/r = 0.900 or 0.752; 
B/p = 33.5 or 57.5. 



(9) 



Figure m shows the result. The agreement [T.N.: be- 
tween the theoretical and numerical relations] is quite 
good, especially when considering that there is no free 
parameter in this comparison, and when noticing that 
the initial state of the clusters, also visible on Figure S] 
was very different from that of the homologous model. 

3.2 Energy distribution 

von Hoerner also gives the energy distribution of the stars 
(v.H., Figure 9 and last column of Table 14). The energy 
is normalized to its mean value E, and transformed into 
another quantity e, for graphical purposes: 



e — arcsinh ( = 
,E 



(10) 



and /(e) de is the fraction of the stars whose energy is 
within e and e -I- de. 

To obtain the corresponding curve in the case of the 
homologous model, one must first calculate the mean en- 
ergy of a star. It is easy to demonstrate that, for a viri- 
alised cluster, this mean energy is minus three times the 
mean kinetic energy (see v.H., Eq. 23), i.e.: 



E = -3-^ = -1.876. 
M, 



(11) 



Furthermore, the fraction of the stars having an en- 
ergy between E and E -f- dE is (Paper I, 2.42): 

FQNiE 

Me ■ 

When using the new variable e, this expression must be 
equal to /(e) de. Thus, using (10), one gets: 



fie) 



FQ' \/E2 + E^ 



(12) 
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Figure 4: Comparison of the density law: homologous model (solid 
line) and numerical experiments {[T.N.: dotted and] dashed lines). 
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Figure 5: Comparison of the energy distribution: homologous 
model (solid line) and numerical experiments {[T.N.: dotted and] 
dashed lines). 



The distribution obtained is plotted in Figure [5l and 
is compared with the resuhs of von Hoerner. Note that, 
here again, there is no free-parameter. The agreement 
is excellent from £ = to e = 2.5, i.e. for the energies 
ranging from to about 6 times the mean energy. This 
interval contains almost all the stars (see v.H., Figure 8). 
In particular, the "shoulder" of the distribution next to 
e = is very well reproduced. For E [sic, e] larger than 
2.5, the results differ; the stars with highly negative ener- 
gies are in excess [T.N.: in the numerical experiments of 
von Hoerner] compared to the homologous model. This 
probably comes from the formation of very close binaries, 
detected by von Hoerner (his Table 5), and predicted by 
theory as described above. 

Let's digress about this topic. One could be surprised 
that the theory leads to the disagreement illustrated in 
Figure [5j although it predicted the existence of binaries. 
In fact, the homologous model presented here, as well as 
all the other theoretical models presented up to now, to 
our knowledge, is fundamentally unable to account for 
binaries. We assume, indeed, that the state of the cluster 
is fully described by the distribution function (Paper I, 
2.1): 

ip{x,y, z,Vx,Vy,Vz,rn,t) or, in short: ip{T,t), 

which gives the probability to have a star in a small vol- 
ume of phase-space, at a given time. Doing this, we im- 
plicitly assume that the stars are independent from each 



other, meaning that the probability of finding a star at a 
given point of phase-space does not depend on the pres- 
ence or absence of stars in other points of phase-space. 
However, this is only a hypothesis and it is not necessar- 
ily true. It seems that this fundamental point has not 
received all the care it deserves; it is even ignored in the 
classical stellar dynamics studies. 

But this hypothesis contradicts the existence of bina- 
ries in the cluster. Indeed, if pairs exist, the presence of 
a star at a point of (ordinary) space leads to a higher 
probability to find another star in the neighbourhood. In 
phase-space, the probability is also modified. To get a 
correct mathematical description of this effect, we must 
introduce the double distribution function: 

(p2{Tl,T2,t), 

which gives the probability that a star is at the point ti 
and another star is at the point T2, at the time t. 

The same way, it would be necessary to introduce a 
triple distribution function ip^ to account for triple stars, 
and so on. These multiple distr ibution functio ns are com- 
mon in plasma physics (see e.g. Delcroixlll963 . Chapter 6) 
and it would be interesting to introduce them in stellar 
dynamics too. This would perhaps allow us to calculate 
exactly, e.g., the formation rate and the properties of bi- 
naries in a cluster. We could also expect from this a more 
precise theory about the effect of encounters. 
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3.3 Evolution 

In Paper I, Chapter VII, we have shown that the ho- 
mologous model is the final state toward which a cluster 
tends, whatever its initial form was. This is clearly con- 
firmed by the results of von Hoerner (v.H., Figures 5 and 
6): his clusters, which initially have a constant density in- 
side a sphere, progressively concentrate in the centre, and 
form a very extended halo at the same time; the density 
law tends toward the theoretical curve of the homologous 
model (Figure S]). 

We have also shown that in a cluster which has ini- 
tially a finite central density, this density rapidly increases 
and tends toward infinity. This phenomenon is clearly 
visible in the Figure 4 of von Hoener, where the "density 
radius" Rp tends to zero. It is easy to verify (v.H., Eq. 10) 
that this radius can become zero only if the central den- 
sity is infinite. One can also directly see the rapid growth 
of the central density (v.H., Table 9). 

Finally, we have seen above that the isolated homol- 
ogous model expands with time (Figure E]). This also 
appears on the right-hand side of the Figure 6 of von 
Hoerner: once the clusters have reached the homologous 
state, the evolution is only a dilatation of the entire clus- 
ter, which is seen on the figure as a translation of all the 
points to the top. This is especially visible in the second 
case. 

Now, let's make a quantitative comparison. We first 
have to adapt the timescales. We have (Paper I, Eq. 8.9): 



dT 



'iM"r'(^)"(5) 



3/2 



(13) 



T being the time variable of the homologous model, while 
t is the real time. We have Ai^ — fnn; and when using 
the definition of the unit of time given by von Hoerner 
(v.H., p. 54), we get: 



— =ni/2 [ln(n)]-i M-'^^(-) 



3/2 



(14) 



The r/R ratio does not vary much with time (v.H., Ta- 
ble 10); we can assume it remains constant and adopt the 
values given by (9). With n = 25, we get: 



dt 
dT 



1.303 or 1.707. 



(15) 



From Paper I, (2.32) and (5.13), we find out: 

ti = -1.188 or - 1.556, (16) 

and (6) gives the theoretical rate of change of the radius, 
at t = 0: 

dlog(0 ^Q 244 or 0.186. (17) 
dt ^ ' 

The observed rate of change for the median radius is 
(v.H., Table 10; we used the mean value, over the entire 
evolution) : 



dlog (rmcd) 
dt 



0.026 or 0.031. 



(18) 



Here, we find a clear disagreement: the two clusters evolve 
much slower than predicted by theory, by a factor of 10 
or 6, approximatively. 

We can also calculate the time required to form an 
infinite central density. Back to the calculation made in 
Paper I, Chapter VII, and assuming that the parameter d 
takes the same value here, we find that a isolated cluster, 
starting from a weakly concentrated state, must get an 
infinite central density after a time: 



I.e. 



T2 = 0.3052, 



t2 = 0.473 or 0.620, 



(19) 



(20) 



from Paper I, (5.12), (2.32) and (15) above. 

But the Figure 4 of von Hoerner shows that the infinite 
central density appears after about t = 6 for the first 
cluster, and t = 4 for the second. Here again, the observed 
evolution rate is smaller than the theoretical one, by a 
factor 13 or 6. 

The reason for such a disagreement has not been 
found. 



3.4 Escape rate 

As mentioned above, the escape rate of an isolated cluster 
is much smaller than that of a non-isola ted cluster, but is 
not exactly zero. It reads (lHenonlll96nl Eq. 15): 



dn 
'dt 



— = -0.00868, 



/ Gmn 



(21) 



where rg is the radius enclosing, in projection, the half of 
the total mass. 

The clusters from von Hoerner have initially the shape 
of an homogeneous sphere of radius unity; in this case, 
we have vq = 0.6083. Furthermore, n = 25. Taking into 
account the choice for the unit of time {[T.N.: v.H.,] p. 
54), we obtain: 



P = -0.0915. 

dt 



(22) 



The total duration of the two examples calculated with 
n = 25 is (v.H., Table 8): 



At ^ 6.23 + 3.86 = 10.09; 



-0.92. 



so that the number of escapers should be An 
In fact, we measure: An = — 1 (v.H., Table 6). 

Four other examples have been calculated with n = 
16. Their duration is not given by von Hoerner. However, 
the total number of relaxation times covered is 88, instead 
of 64 for the examples with n = 25. The number of es- 
capers should be larger in proportion, i.e. Art = —1.27. 
The number measured is: An = — 1. We note that here, 
the agreement between the theory and the numerical ex- 
periments is very good. 
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4 CONCLUSION 

We have compared the results obtained for the same prob- 
lem: structure and evolution of an isolated cluster, thanks 
to two very different methods, almost opposite: the first 
one (the homologous model) describes the cluster with a 
continuous distribution function, whose evolution is ruled 
by a Fokker-Planck type equation; the second one (nu- 
merical experiments) represents the cluster as a discrete 
system made of material points, whose evolution is given 
by the ordinary laws of mechanics. 

The results are generally in very good agreement. This 
constitutes a solid argument in favor of the validity of 
both theoretical methods. 

Some hypotheses are also confirmed. In the case of the 
homologous model, we had to assume an isotropic distri- 
bution for the velocities everywhere in space, in order to 
reduce the complexity of the equations. No restriction 
of this kind exists in the numerical experiments, where 
the entire evolution of the system of points is computed 
exactly. The agreement of the results, particularly for 
the density law (Figure 2]), indicates that this hypothesis 
about isotropy does not strongly affect the structure of 
the cluster. We can reasonably suppose that the same 
applies to the non-isolated clusters, studied in Paper I. 

Furthermore, in the case of the numerical experiments, 
one could ask whether a system of only 25 points can 
properly represent a real cluster, made of a much larger 
number of stars. But the homologous model, with its con- 
tinuous distribution function, can be considered as the 
limiting case of a system made of an infinite number of 
points. The similar results obtained in both cases are 
thus, very probably, also valid for all the intermediate 
cases, i.e. from n — 25 to n — od; this range includes all 
the real stellar systems, which are thus "surrounded" by 
the two methods. We can also claim that, for the ques- 
tions we address, 25 is a good approximation of infinity. 

Therefore, it seems recommended to develop further 
the method of the numerical experiments, which is still 
in its early stages. Addressing simultaneously the ques- 
tions of stellar dynamics with theory and numerical ex- 
periments would, certainly, allow us to make more solid 
and faster progresses than in the past. 
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